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We present a combination of density functional theory with dynamical mean field theory to study 
correlation effects in organic layered charge-transfer salts. For that purpose, we propose a scheme to 
construct molecular Wannier functions using atomic orbitals as starting point. We compute spectral 
and optical properties of k-(BEDT-TTF)2Cu[N(CN)2]C1 and demonstrate from first principles that 
interdimer transitions are strongly affected by correlations, leading to a Hubbard peak in the optical 
conductivity, while intradimer transitions are uncorrelated, giving rise to a peak consistent with LDA 
predictions. We discuss our results in the context of reported infrared reflection measurements. 
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Electron correlated materials show a rich phase dia- 
gram originating from the interplay of factors like mag- 
netic ordering, frustration, doping, pressure, or temper- 
ature. Organic charge transfer (CT) salts are unique 
correlated systems with a complex phase diagram that 
can be finely tuned by chemical substitution or moder- 
ate pressures. Observed phases in these systems include 
metals, Fermi liquids, Mott insulators, antifcrromagnets, 
spin liquids, and unconventional superconductors!^ 

An extensively studied class of CT salts are based 
on the molecules bis-(ethylenedithio)tetrathiafulvalene 
(BEDT-TTF, or shorter ET) crystallized in the k phase. 
In these k-(ET)2X salts, electron donors (ET) and elec- 
tron acceptors (X) form alternating layers, with two ET 
molecules facing each other in ET dimers; these dimers 
are arranged in a triangular lattice. For monovalent an- 
ions X one electron is transferred from each dimer (ET)2 
to each anion formula unit. In effect, the dimer layers 
possess one hole per dimer and their highest occupied 
molecular orbitals (HOMO) give rise to a half-filled con- 
duction band. Hence, band structure calculation^E pre- 
dict the dimer layers to be metallic. However, the exper- 
imentally observed ground state depends on the choice 
of the anion: even for the example of the isostructural 
compounds k-(ET) 2 Cu[N(CN) 2 ]C1 (in short K-Cl) and k- 
(ET) 2 Cu[N(CN) 2 ]Br (k-Bi), the ground state can be as 
different as a Mott insulator for k-C1 and a Fermi liquid 
for K-Br at low temperatures and ambient pressure!^ In 
fact, k-CI can be driven through the insulator-to-metal 
transition (MIT) by the gradual substitution of CI for 
isovalent Br which changes the volume of the system and 
thereby decreases the ratio U/W of the Hubbard U to 
the band width W . Being close to the MIT, K-Br is a 
'good' (coherent) metal with a Drude peak only below a 
coherence temperature T* sa 50 K. At higher tempera- 
tures k-C1 is a semiconductor with a gap of E g = 800 K 7 
while K-Br shows 'bad metal' behavior with strong scat- 
tering preventing coherent transport and suppressing the 
Drude peak. Theoretically, this crossover from a Fermi 
liquid at low T to the 'bad' metal state at higher tem- 
peratures in K-Br is obtained from dynamical mean field 



theory (DMFT) calculations where gradual destruction 
of coherent quasiparticles occurs with increasing temper- 
ature,^ albeit at higher T* than observed experimentally. 

The fact that small chemical modifications lead to 
qualitative changes in behavior together with the impor- 
tance of electronic correlations in these materials make 
it clear that a realistic description requires both (i) de- 
tails of the band structure as well as (ii) a proper treat- 
ment of strong correlations. However, many-body stud- 
ies of the k-(ET) 2 X salts have so far been limited to 
minimal model calculations^^! f the Hubbard or ex- 
tended Hubbard Hamiltonian on an anisotropic triangu- 
lar latticed In this work, we go beyond model calcu- 
lations and present a combination of ab-initio density 
functional theory calculations in the full potential lin- 
earized augmented plane wave (FLAPW) framework^ 
combined with DMFT^H (LDA+DMFT). With this ap- 
proach, the above requirements for a realistic description 
of these systems are fulfilled to a large extent. To our 
knowledge, this is the first LDA+DMFT calculation on 
an organic molecular crystal. Specifically, we focus on 
the spectral and optical properties of k-C1 at room tem- 
perature. While the kinetic part of the Hamiltonian is 
described with full DFT precision, we use a few values of 
the Hubbard interaction strength U in order to illuminate 
the effects of different U/W ratios. We will show that, 
whereas it is well known that the DMFT method does 
not yield a perfect description of the underlying Hubbard 
Hamiltonian on a triangular lattice in two dimensions, 
this approach provides insight into the contributions to 
the optical conductivity at an unprecedented level. Our 
calculations demonstrate that interdimer and intradimer 
transitions are responsible for two principal features ob- 
served in optical conductivity measurements at low tem- 
peratures and identify the interdimer feature to be re- 
lated to correlations. These results corroborate previ- 
ous conjectures based on phenomenological and minimal 
model assumpt ions ! 1 3 * 1 ^ 

In order to perform LDA+DMFT calculations on 
molecular crystals we have to construct localized Wannier 
functions for molecular orbitals as required by DMFT. 
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The construction method based on projecting Bloch 
states onto pur e atomi c orbitals with subsequent or- 
thonormalizatior l 12 * 15 * 16 ! is fast and stable. In its stan- 
dard form, though, the projection method is, by defini- 
tion, not designed for molecular orbitals. In this work, we 
propose a scheme to construct molecular Wannier func- 
tions using atomic orbitals as a starting point. The key 
element of this scheme is the diagonalization of the oc- 
cupation matrix written in the basis of atomic orbitals 
within the subspace of correlated bandsP^The real space 
representation of the resulting dimer HOMO Wannier 
function of re-Cl (based on the crystal structure reported 
in Ref. 1181 with space group Pnma) is shown in Fig. [I] 




FIG. 1. (Color online) Structure of k-(ET) 2 Cu[N(CN) 2 ]C1 
seen along the c axis with a Wannier function corresponding 
to the bands crossing the Fermi level. 

For solving the impurity problem we employed the hy- 
bridization expansion continuous-time quantum Monte 
Carlo method^ as implemented in the ALPS cod e 1 20 1 21 1 
where we used 2x 10 7 Monte Carlo sweeps throughout our 
calculations at an inverse temperature (3 = 40 cV _1 , cor- 
responding to room temperature. Note that in the crys- 
tallographic unit cell, there is a manifold of four bands 
around Ep as the crystallographic unit cell contains two 
organic layers with a total of four dimers due to the pres- 
ence of the anion. The dimers are equivalent, though, 
so that only one dimer HOMO has to be considered in 
our single-site DMFT; the other three orbitals are re- 
lated by symmetry. Since the quantum Monte Carlo al- 
gorithm operates on the imaginary frequency axis, the 
calculation of dynamical quantities like spectral func- 
tions and optical conductivity requires analytic contin- 
uation to the real axis. We performed stochastic ana- 
lytic continuation^! on the self energy for obtaining the 
spectral functions and directly on the optical conduc- 
tivity <j{iv) for the calculation of optical properties. 23 
Concerning appropriate interaction parameters for these 
systems, considerably distinct values for the Hubbard U 
have been reported. In this work we analyze the effect 
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FIG. 2. (Color online) Momentum-resolved spectral function 
for U = 0.6 eV. 



of the values U = 0.5/0.6/0.7/0.84 eV. This range is 
guided by two estimates: U — 0.85 eV obtained for a 
similar but arguably more strongly correlated compound 
k-(ET) 2 Cu 2 (CN)3 from constrained random phase ap- 
proximations', and U « 0.27 eV extracted from optical 
conductivity measurements^ and model considerations.^ 

In Fig. [2] we show the calculated band structure of 
k-CI in form of the momentum-resolved spectral func- 
tion for U = 0.6 eV along with the LDA band energies. 
The bands at the Fermi energy E-p — originate from 
the interdimer hopping, in particular hopping between 
dimers on the same layer; the interlayer hopping is very 
small, so that the four bands are composed of two al- 
most degenerate pairs of bands. The correlation in band 
space acts almost exclusively on these bands, splitting 
them into renormalized excitations of (mass-enhanced) 
quasiparticles and a spectral weight transfer to an up- 
per and lower Hubbard band which manifests itself as 
blurry background. On the other side, the charge trans- 
fer between the ET molecules within a dimer (intradimer) 
induces the splitting between the bands right below the 
correlated manifold and the correlated bands. For uncor- 
rected states, the excitations basically lie on top of the 
respective LDA energies which validates the accuracy of 
our analytic continuation. 

The development of the Hubbard bands with U is vis- 
ible from the momentum-integrated spectral functions in 
Fig. [3| At Ep , the spectral functions are suppressed com- 
pared to LDA and Hubbard bands are formed for all val- 
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ues of U, getting more pronounced as U is increased; the 
maxima of the Hubbard peaks are located at uj ~ ±U/2, 
as expected. However, a quasiparticle peak at Ep re- 
mains present up to U — 0.7 eV; the largest value of U 
shown, U = 0.84 eV, is just on the verge of opening a gap 
at Ep. It is known, however, that the critical value of U 
for the MIT, U c , is significantly overestimated by DMFT 
in low dimensions: for the triangular latice in two dimen- 
sions as m our case, DMFT yields U c « 15 \t^ whereas 
extrapolated cluster approximations^ yield U c ~ 7\t\. 




FIG. 3. (Color online) LDA density of states and the 
LDA+DMFT spectral functions for different values of U. 

Next, we investigate the optic al properties of k-C1 
where, like in previous works JEliiJ we concentrate on the 
polarization E || c, i.e., along the linear chains in the 
triangular lattice. The authors of Ref. [13] performed a 
comprehensive study on the evolution of the optical con- 
ductivity spectra from infrared reflection measurements 
on K-(ET) 2 Cu[N(CN) 2 ]Cli_ x Br x tuned through the MIT 
by increasing the Br concentration x. At room tempera- 
ture, they observed a broad midinfrarcd absorption peak 
between 1600 cm -1 and 4200 cm -1 , in agreement with 
previous optical studies on k-C1 and K-Br (see Ref. 1 131 and 
references therein). At low temperature, a Drude peak 
evolves for the compounds with high Br concentration 
which marks the onset of metallicity at x ~ 0.7 whereas 
no Drude peak is visible for lower Br content indicating 
an insulating state without coherent quasiparticles. 

Importantly, at low temperatures the broad midin- 
frared peak splits into two peaks in the pure CI and low 
Br concentration compounds, fitted by two Lorentzians 
at «2200 cm' 1 and 3200 cm -1 ; for high Br content this 
splitting is very weak. From this doping dependence, 
it was concluded^ that the first peak is a correlation- 
induced feature due to electron transitions between the 
lower and upper Hubbard bands, while the second 'dimer 
peak' was assigned to the intradimer charge transfer. 
For the metallic Br rich compounds, the spectral weight 
transfer from the Drude peak to this Hubbard peak is 
small and the Hubbard peak is suppressed or less pro- 
nounced. Since the Hubbard peak is expected to ap- 



pear at uj ss U, this would determine the effective U as 
U « 2200 cm" 1 = 0.27 eV. 
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FIG. 4. (Color online) Calculated LDA+DMFT optical con- 
ductivity for U = 0.6 eV together with the LDA result and 
the experimental curve at T = 300 K. 

Our method allows us to test this interpretation; as the 
energy window of our LDA+DMFT calculation contains 
both the correlated manifold at the Fermi energy as well 
as uncorrelated bands away from Ep, all transitions can 
be inspected on equal footing. As mentioned above, the 
symmetry of the organic layers is higher than the sym- 
metry of the crystal, therefore the four bands around Ep 
can be considered to originate from one band which is 
folded due to the unit cell that is twice as large as re- 
quired by ET layer symmetry. These bands thus create 
the intraband contribution to the optical conductivity 
discussed in the following; physically they correspond to 
the mierdimer hopping. All other transitions are termed 
mferband transitions; among others, these contain tran- 
sitions related to the miradimer charge transfer. 

The LDA+DMFT calculated optical conductivity for 
U = 0.6 eV together with the experimentally measured 
spectrum of K-Cl at T = 300 K is shown in Fig. [4] 
Overall, the calculated total optical conductivity a to t fea- 
tures one dominant peak at approximately 3450 cm' 1 , 
i. e., close to the experimental peak position (we do not 
consider vibronic modes) and a second much higher- 
frequency peak between 8000 cm' 1 and 17000 cm' 1 . 
The peak at higher frequencies originates from interband 
transitions between uncorrelated states and hence corre- 
sponds in frequency with LDA (also shown in the figure); 
the analytic continuation smears out the fine structure 
of this peak. The position of the low-frequency peak of 
atot is roughly centered at the same position as the LDA 
results, but it is strongly enhanced in spectral weight, 
in accordance with experiment. While our calculations 
are at room temperature, they already capture the ef- 
fects of correlation observed more pronounced in mea- 
surements at low T. Decomposed into intra- and inter- 
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band contributions, we find that erintcrband roughly co- 
incides with the LDA results in position and spectral 
weight. In addition, we find an intraband absorption 
at w«3550cm _1 «0.75 U which we identify as intraband 
Hubbard transitions. While one expects a peak at uRiU 
from transitions between the Hubbard bands, a U/2 peak 
mixes in since a quasiparticle peak is still present in our 
calculations (for the chosen U = 0.6 eV) and allows for 
transitions between the quasiparticle peak and the Hub- 
bard bands. Accordingly, our data show a remnant of 
the Drude peak not observed experimentally for the pure 
CI compound. The previously mentioned limitations of 
DMFT here force us in our choice of U values to compro- 
mise between the correct semiconducting phase (favored 
by a large U) and the correct Hubbard peak position 
(favored by a moderate U). 




FIG. 5. (Color online) Optical conductivities for different 
values of U . The peak positions of the intraband (interband) 
contribution are: u) » 3000 (3200) cm" 1 for U = 0.5 eV, ui « 
3550(3200) cm -1 for (7 = 0.6 eV, and ui » 4250 (3500) cm -1 
for U = 0.7 eV. The arrows show the respective peak positions. 

As for the nature of the absorption peaks, Fig. [5] 
shows their evolution with U. The intraband contribu- 
tion to the conductivity moves in frequency proportional 
to U, with the peak position consistently corresponding 
to ss 0.75/7 while the interband absorption is largely in- 
sensitive to U . This analysis demonstrates the correlated 
nature of contraband and uncorrelated nature of cr interband . 
In all cases, the two peaks are too close to each other to 
be individually resolved from the total optical conduc- 
tivity alone. Furthermore, Fig. [5] shows how the Drude 
peak develops for decreasing U values; Fig. [5] (a) with 
U = 0.5 eV is actually similar to the case of K-Br as 
analyzed in Fig. 8 of Ref. [131 

Quantitatively, the suppression of the Drude peak as 
a function of U and the redistribution of the intra- 
band spectral weight is presented in Fig. [6] where we 
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FIG. 6. (Color online) Distribution of the intraband spectral 
weight in terms of the effective number of charge carriers. 



follow Ref. [5) and plot the integrated spectral weight 
Jo contraband (a/ )cfc/ representing the effective number of 
charge carriers N e g. In this representation, the number of 
charge carriers in LDA by definition equals the number of 
conduction electrons, i. e. one, and all the weight is con- 
centrated in the (infinitesimally narrow) coherent Drude 
peak which is only broadened by temperature. Upon in- 
clusion of correlations, the kinetic energy of the electrons 
is diminished, which corresponds to mass enhancement 
(in Fermi liquid theory) or a reduction of the number of 
effective charge carriers as we observe. 

Summarizing the analysis of the optical conductiv- 
ity calculations, we find agreement with the princi- 
pal assessment of Ref. 13, namely that the two finite- 
frequency peaks in the experimental optical conductiv- 
ity of k-C1 at low temperatures originate, respectively, 
from correlation-induced intraband (intcrdimer) contri- 
butions ('Hubbard peak'), scaling with U, and interband 
(intradimer) transitions which are unaffected by correla- 
tions ('dimer peak'). The fact that DMFT overestimates 
the critical U of the triangular lattice in two dimensions 
forces us to describe k-C1 with a somewhat high interac- 
tion strength of U=0.6 eV, suggesting that future stud- 
ies combining LDA with many-body techniques beyond 
DMFT like cluster extensions of DMFT^l may further 
improve the quantitative agreement with experimental 
measurements. This is beyond the scope of the present 
study. 

In conclusion, we presented the first LDA+DMFT 
study on the spectral and optical properties of the organic 
CT salt k-(ET) 2 Cu[N(CN) 2 ]C1. Our results provide 
an ah initio-h&sed theoretical evidence for the double- 
natured origin of the infrared peak in the optical con- 
ductivity of k-(ET) 2 Cu[N(CN) 2 ]C1. We could identify 
intraband transitions within the correlated manifold and 
interband transitions due to charge-transfer within an ET 
dimer. The proposed projection method for construct- 
ing non-atom-centered Wannier functions in the FLAPW 
framework is computationally efficient and can be applied 
to a great variety of correlated organic as well as inor- 
ganic systems with (quasi-)molecular orbitals. 
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